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ABSTRACT 


Computational predictions and experimental measurements were 
made on the Naval Postgraduate School’s cold-flow turbine test rig. The 
test turbine was the Space-Shuttle Main Engine, High-Pressure Fuel 
Turbopump, Alternate Development Model, designed and manufactured by 
Pratt & Whitney. The flow-field around the first-stage rotor end-wall 
region was measured using a laser-Doppler velocimetry (LDV) system. 
Measurements were taken at two axial locations over the rotor blade tip 
and at three radial locations from the end-wall casing. Three 
circumferential velocity profile measurements were taken downstream of 
the first-stage using a three-hole pressure probe. All measurements were 
taken at a referred rotational speed between 4781 and 4904 rpm. A 
computational fluid dynamics model of the combined first-stage stator and 
rotor was developed. Predicted velocity data from this model were 


extracted for comparison to the rotor exit plane probe measurements. 
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I. INTRODUCTION 


A. PURPOSE 


This thesis describes experimental flow field measurements in the 
first-stage turbine rotor end-wall of the Space Shuttle’s Main Engine 
(SSME) High-Pressure Fuel TurboPump (HPFTP). It also describes the 
generation of a combined turbine stator and rotor computational fluid 
dynamics model. Comparisons were made between the experimental 
measurements and the numerical predictions of the rotor exit plane flow 


field. 


B. OVERVIEW 


With the advancement of computers in the past decade, advanced 
viscous flow fields can be computed within a reasonable period of time. 
The use of computational fluid dynamic (CFD) codes can save money and 
time in the development of turbomachines. With the help of these codes 
more efficient blade shapes can be designed to decrease aerodynamic 
losses and heating. Such designs can result in longer service life and 
lower operational costs of the gas turbine engine. 

The Naval Postgraduate School’s Turbopropulsion Laboratory (TPL) 
has been supplied with the Alternate Turbopump Development (ATD) 
model of the SSME HPFTP, designed by Pratt & Whitney. In Figure l,a 
schematic of the SSME from Sutton [Ref. 1: 1992], the HPFTP is depicted 
on the left-hand side of the diagram. This device consists of a two-stage, 
axial flow turbine. The turbine is driven by the flow of hydrogen gas and 
steam, which is the product of the combustion of hydrogen and oxygen. 
The turbine in turn drives a 3-stage hydrogen pump. The HPFTP was 
designed to operate at an inlet temperature and pressure of 1900 degrees 


Rankine (°R) and 5200 pounds per square inch (psia), while developing 


73,000 horsepower (hp) at a rotational speed of 40,000 revolutions per 


minute (rpm). 
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Figure 1. The Space-Shuttle Main Engine 


The HPFTP has been modified for the Turbine Test Rig (TTR) to 
facilitate LDV measurements for the purpose of validating viscous flow 
codes. Several students have undertaken this project. The first was 
Studevan [Ref. 2: 1993] who initially designed and installed the cold-flow 
test facility. Rutkowski [Ref. 3: 1994] redesigned parts of the bearing 
housing and was responsible for the initial numerical modeling of the 


turbine stator. Greco [Ref. 4: 1995] further modified the TTR, adding the 
2 


data acquisition system, and continued the numerical analysis of the first 
Stage stator and rotor. Southward [Ref. 5: 1998] redesigned the layout of 
the TTR to accommodate the LDV system and performed the first LDV 


measurements in the turbine. 


1. Experiment 


The ultimate goal of this research was to collect LDV measurements 
to validate Computational Fluid Dynamics computer models of the cold- 
flow immediately adjacent to the end-wall of the turbine in the first-stage 
rotor tip clearance region. LDV measurements were taken at three axial 
positions with three radial depths for each axial position. Cobra probe 
surveys at three circumferential positions were also taken at the first- 


stage rotor exit plane. 


Z Numerical Simulation 


Three different numerical simulations were conducted. The flow 
field of the first-stage stator was investigated using RVC3D (Rotor 
Viscous Code 3-D) [Ref. 6: 1992], and using SWIFT [Ref. 7: 1997]. The 
combined first-stage stator and rotor flow fields were investigated using 
the SWIFT code. Two grid generation programs were used. The Stator 
and rotor grids were first developed in TCGRID (Turbomachinery C-Grid) 
[Ref. 8: 1990] and later developed in TCGRID version 202 [Ref. 9: 1996], 
which is compatible with the SWIFT code. 


Ja Past and Ongoing Research 


Several studies on computational studies of flow through turbines 
and LDV measurements in tip clearance areas have been conducted over 
the past few years. Most recently Ameri, Steinthorsson and Rigby [Ref. 
10: 1998] used computational models to study the effect of tip clearance 
and casing recess on heat transfer and stage efficiency in axial turbines. 


However they only computed a rotor flow field without an upstream stator 


2 


and presented their numerical predictions without comparison to any 
experimental data. Stahlecker and Gyarmathy [Ref. 11: 1998] conducted 
3-D LDV measurements while investigating turbulent flow in a centrifugal 
compressor vaned diffuser. This study is the most recent in the ongoing 
effort to acquire three-dimensional LDV measurements in turbomachines. 
For a more complete list of references pertaining to turbomachinery LDV 
measurements, particularly those in turbines, refer to Southward [Ref. 5: 


LOO e 


Il, EXPERIMENTAL SETUP 


A. TEST FACILITY AND HIGH-PRESSURE FUEL TURBOPUMP 
TURBINE 


The Naval Postgraduate School’s High Speed Turbopropulsion 
Laboratory (TPL) was the site of all testing done on the High Pressure 
Fuel Turbopump (HPFTP). The cold-flow Turbine Test Rig (TTR) was 
driven by compressed air supplied by an Allis-Chambers compressor. A 
schematic of the air supply system 1s shown tin Figure 2 and the 


compressor is shown in Figure 3. 
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Figure 2. Schematic of Allis-Chalmers 12 Stage Axial Compressor 





Figure 3. 12 Stage Axial Compressor 


Compressed air was routed to the test cell through large pipes. The flow 
of compressed air to the test cell was controlled using one in-line and two 
dump valves, which were adjusted manually at a console outside the test 
cell (Figure 4). The test cell itself housed the Turbine Test Rig (TTR), 
laser-Doppler velocimetry (LDV) system and Scanivalve. A schematic of 
the TTR is shown in Figure 5 and depicted in Figure 6. The bearing 
housing was rebuilt after bearing failure and the throttle guide rod, 
movable back pressure plate and shaft cover were removed for ease of 


maintenance. 
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Figure 4. Photo of Dump-valve Controlls for Compressed Air Control 


The data acquisition system was housed in the main room of the 
TPL. Two personal computers (PCs) were used, one to collected data 
from the TTR and the other collected data from the LDV system. Figure 7 
‘shows the layout of the test cell. For a detail description of all 


experimental equipment refer to Southward [Ref. 5: 1998]. 
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Figure 6. Photo of Turbine Test Rig : 
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Turbine Test Rig Test Cell, from Ref. [5] 


5: COBRA PROBE SETUP 


A three-hole (Cobra) pressure probe was used to measure total 
pressure, Mach number and flow angle at three circumferential positions 
of the first-stage rotor exit plane. The Cobra probe’s radial and yawing 
movements were controlled from the main room of the TPL using a 
control box which housed a DC-to-AC power converter and the actuators 
for the Motion system. The three-hole probe was inserted through an 
access hole and attached to the TTR using a mounting plate. The access 
hole was adjusted to the three survey sites by detaching the aft portion of 
the TTR and rotating it to the desired position and then reattaching. 


Figure 8 shows the Cobra probe attached to its calibration stand. 





Figure 8. Cobra Probe Mounted on Calibration Stand 


Two pressure transducers were connected to the probe. The center hole 


was connected to a single transducer and used for total pressure 
10 


measurement. The two outboard ports were connected to either side of a 
single differential pressure transducer. The Cobra probe and associated 
pressure ports were connected to a Scanivalve by plastic tubing. This 
tubing was replaced with more ridge plastic tubing to correct the pressure 


lag discovered in previous experiments. 


C. LASER DOPPLER VELOCIMETRY SETUP 


The LDV system was set up to measure the axial and 
circumferential components of flow in the TTR’s first-stage rotor end- 
wall region. An argon-ion laser and associated optics were mounted on a 
traverse table system (Figure 9). The green and blue beams produced by 
the laser were split and one side of each beam was sent through a 
frequency shifter before being focused inside the TTR observation port. 
The frequency shifters were found to be faulty and sent out to be 


recalibrated. 





Figure 9. Laser and Optics on Traverse Table 
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Two six-jet atomizers were used to seed the flow with anhydrous glycerin. 
The seeding material was introduced into the flow with a modified wand 
far upstream and adjacent to the wall (Figure 10). The seeding wand was 
redesigned to introduce the seeding particles closer to the turbine wall for 
better seeding density. The interference of the laser beams produced 
fringes. The reflected light, produce by the seeding particles passing 
through the fringes, was collected through receiving optics by two 
photomultipliers. The signals from the photomultipliers were sent to an 
IFA750 for correlation with the once per rev output of the rotating 
machinery resolver. A TSI parallel interface card, mounted in a 386 
personal computer (PC), collected the information from the IFA750. 
Phase Resolved Software version 2.06 was used to analyze the data 


collected from the IFA7S0O. 
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Figure 10. Two Six-Jet Atomizers 
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D. TTR DATA ACQUISITION 


The TTR data acquisition system was located outside of the test cell 
in the lower portion of the control room (Figure 11) and consisted of a 
486 PC with interface boards, Signal Conditioner Unit (SCU), and three 
Hewlett-Packard data acquisition modules. The Scanivalve was located in 
the test cell (Figure 12). Data acquisition was controlled by a LabView 
Software Package. For a detailed description of each component refer to 


Southward [Ref 5: 1998]. 





Figure 12. Scanivalve 
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KE. EXPERIMENTAL PROCEDURE 


Two different experiments were conducted on the TTR. After the 
replacement of the bearings in the bearing housing, Cobra probe surveys 
were conducted at three circumferential locations in the rotor exit plane. 
LDV measurements were then conducted in the end-wall region of the 


first-stage rotor. 


ie Cobra Probe Measurements 


A three-hole pressure probe was used to measure the exit flow 
conditions of the first-stage rotor. Exit total pressure, flow angle and 
Mach number were measured. The probe was located .31 inches aft of the 
rotor tip trailing edge. Four complete surveys were conducted from the 
end-wall to the hub at three circumferential positions, with two taken at 


position 2. Figure 13 shows the circumferential position of the surveys. 
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Figure 13. Positions of Cobra Probe Survey 
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Set up of the three-hole pressure probe is described in Southward [Ref. 5: 
1998, pgs 44-45]. Probe depth was monitored by the LabView program 
ACTUATOR.VI and controlled by a switch box which was manually 
operated from the control room’s data-acquisition station. After setting 
the probe to the desired depth the flow angle was determined by nulling 
the signal from the differential pressure transducer. The probe was 
rotated left or right until the out-board pressures were equal. This was 
indicated by a near zero read-out from the digital volt meter (DVM). 

After nulling was accomplished, the data acquisition program 
VEL_PRFL.VI was run to collect pressure and temperature readings from 
the various probes on the TTR. The program then calculated and 
displayed Mach number. The VEL_PRFL.VI program produced an output 
file named VEL_PRFL.DAT. The data from this file were extracted and 
plotted with the Microsoft Excel program. Probe surveys were first 
conducted at position 2 (Figure 13). The probe was then removed and the 
aft casing rotated to position 3. The pressure probe was then reinstalled. 
Each time the probe was removed, care had to be taken in order not to 
scrape the probe against the walls of the access port. After completion of 
the survey at position 3, the probe was moved to positions 2 and 1 for 


surveys at each. The probe was then removed from the TTR. 


ee LDV Measurements 


All LDV measurements were taken close to the first stage rotor end- 
wall. Southward [Ref. 5: 1998, pgs 41-54] reported procedures for start- 
up, shut-down, and positioning of the LDV systems. LDV surveys were 


taken at the first two axial locations (Table 1). 






Distance from center ais 
LDV hole (inches) 


Table I. Axial Position of LDV holes 





Three radial depths were surveyed at each axial location. Table 2 lists the 
depths, in inches from the end-wall and percent span of the rotor, for each 


radial location. 


Distance from 
Radial Position Outer casing end-wall 
(% Span of rotor) 





Table 2. LDV Radial Positions 


After the laser was moved to the desired location for measurements, the 
circumferential laser beam was shifted by 10MHz to eliminate blade 
backscatter and the outputs from both photomultipliers were observed 
using the PHASE resolve software running on a PC. Appendix I of 
Southward [Ref. 5: 1998, pgs 131-132] lists the menu settings for the 
PHASE resolve software. All settings were identical with the exception 
of the Automatic Filter Selection, which was 5-30 MHz for both channels. 
The PHASE software enabled the data to be taken every 0.1°. First-stage 
rotor blades were spaced every 7.2°, equating to 50 blades. Data were 
taken every 0.1° around the first-stage rotor at the forward axial position. 
The data acquisition window for the center axial position was adjusted to 


blank out passage of the blade through the probe volume. Again, data 
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were taken around the full 360° of the rotor, but only at positions between 
the blades. 

Data were reduced using the PHASE resolve software. All the bins 
of data were combined into one blade passage. For the forward axial 
position the bins were averaged into 72 bins, corresponding to one bin 
every 0.1°. The center axial position was window-averaged into 30 and 
18 bins, equating to 3.0° and 1.8° respectively. The output files from the 
data reduction were then converted to a usable format using the 
‘“Phase3.for” Fortran program modified by Southward [Ref. 5: 1998]. 


Data were then plotted using the Microsoft Excel program. 


we 





Il. COMPUTATIONAL FLUID DYNAMICS 


A computational fluid dynamics (CFD) model was constructed to 
simulate the flow through the first-stage of the SSME HPFTP ATD for the 
purpose of comparing experimental results with the numerical model. The 
geometry of the first stage stator and rotor were obtained by Greco [Ref. 
4: 1995]. Greco produced the first-stage stator and rotor grids using 
Turbomachinery C-Grid (TCGRID), [Ref. 7:1990]. The flow was then 
simulated using Rotor Viscous Code 3-D (RVC3D), [Ref. 6: 1992]. The 
stator and rotor models were run separately through RVC3D. 

The purpose of this study was to run the stator and rotor together 
using a new program under development by Dr. Chima of NASA Lewis 
Research Center (LERC) called SWIFT, [Ref. 9: 1997]. This new 
program allows for multiple grids to be combined and rotated with respect 


of one another. 


A. MODIFIED STATOR GRID GENERATION 


The old version of TCGRID was used to develop the modified stator 
grid. Since there was no surviving data on disk of Greco’s original stator, 
the data were scanned and an optical character recognition (OCR) program 
was used to convert the scanned image into editable text. This text was 
then compared to the original and errors in the text corrected. Greco’s 
stator grid was then modified to model a step in the end-wall casing just 
aft of the stator (Figure 14). The height of the step was represented by 
ten grid points and this distance was represented by multiple radii to 
curve the edges of the step. The new wall was then edited into the 
“stator.in” file for input into TCGRID. The output from TCGRID was 
then run with RVC3D and SWIFT to compare the results of the two CFD 


programs. Appendix A shows the results of these runs. 
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Figure 14. Step in End-Wall 


B. COMBINED GRID GENERATION 
The combined grid generation used TCGRID version 107 [Ref. 9: 


1997]. Because the data received by Greco were for the hot-flow, and the 
TTR is a cold-flow turbine, he applied a thermal! shrink factor of 99% to 
all radial and chord-wise dimensions. Five radial blade surfaces, along 


with modified hub and end-wall surfaces, were included in the “stator.in” 


and “rotor.in” files for TCGRID. 


ie Stator Grid Generation 


The first-stage stator grid used the same five blade surfaces as 
Greco’s stator. The hub and end-wall were modified to account for a new 
position of the stator exit plane and the modeling of a step located in the 
rotor end-wall. Both the stator and rotor needed full hub and end-wall 
geometry inputs in order to match the k-planes of their grids. The stator 
inlet plane was at an axial location of zero. The exit plane was moved to 
an axial position of 1.5 inches so that the step geometry would be 
included in the rotor grid. The coordinate system used for a general fitted 


body is shown below in Figure 15. 
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Figure 15. General Fitted Body System, from Ref. [6] 


TCGRID version 107 produced grids in the PLOT-3D format and also 
added a dummy grid line in the j-direction so that the grids could over lap 
by one cell. The stator grid input dimensions were 135x31x57. TCGRID 
produced a grid of 135x32x57. The difference in the j-direction accounts 
for the added dummy grid line. The tip gap grid clustering of the rotor 
had to be included in the “stator.in” input file in order to match the 
Stator’s exit k-planes with the inlet k-planes of the rotor. FAST was use 
to view the output of TCGRID. Figures 16 and 17 show the generated 


Stator grid. 
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Figure 16. 3-D First Stage Stator Grid (135x32x57) 
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Figure 17. 2-D Stator Hub C-grid (k=1) 
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ma Rotor and Tip Grid Generation 


The first-stage rotor grid was produced using TCGRID version 107. 
TCGRID produced the rotor C-grid and the tip O-grid together. A total 
grid of size 235x31x57 was specified. The produced grid was of size 
235x32x57. Again, the difference in the j-direction accounted for the 
added dummy grid line. A tip gap of 0.045 inches was also specified 
which gave the tip O-grid a dimension of 113x13x13. The rotor accounted 
for k-dimension of 1-45 with the tip filling in from 45-57. The rotor 
inlet-plane was located at an axial position of 1.5 inches. The rotor’s 
exit-plane was located at the probe survey location of 2.65 inches. 


Figures 18 through 20 show the generated rotor and tip grids. 
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Figure 18. 3-D First-Stage Rotor Grid (235x32x57) 
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Figure 19. 2-D First Stage Rotor Hub (k=1) 
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Figure 20. 3-D Tip O-grid (113x13x13) 
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ae Stage Grid Generation 


After the stator and rotor grids were complete, several checks had 
to be made to ensure that the grids were compatible. The rotor grid had to 
overlay the aft portion of the stator grid by exactly one cell. To this end 
the rotor inlet was squared by setting ‘rcorn’ equal to zero (in the input 
file to TCGRID). The one cell overlap was met by setting ‘dsmax’ (the 
furthest most grid spacing away from the blade) of the rotor equal to 
‘dswte’ (the spacing in the wake at the exit) of the stator. The grids were 
then viewed in FAST to ensure that all the k-plans matched and that there 
was exactly one cell overlapping. The grids were combined using a 
program called MULTIX to combine the grids. Figure 21 shows the first- 
stage stator and rotor. The combined grid output from this file was 


renamed to “fort.1” to be read by the flow solver SWIFT. 
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Figure 21 Combined First-Stage Stator and Rotor 


The program MGRID was written to read the single multigrid “fort.1” file 
and reproduce it three times. This offered a better view of the first-stage 


stator and rotor and is shown in Figure 22. The one cell overlap can be 
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seen Clearly in this figure. Appendix B contains the Fortran code for the 


MGRID program. 
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Figure 22. Combined 3-D First-Stage Stator and Rotor Grids 


C. COMPUTATIONAL SCHEMES 


Two flow solvers were used. RVC3D was first run on the modified 
stator grid describe previously. SWIFT was then tested on the modified 
stator. After verification by RVC3D that the code was running properly, 
SWIFT was used to model the flow through the stage. The Naval 
Postgraduate School’s Cray Y-MP EL-98 supercomputer was used for all 


computations. 


1. RVC3D 


RVC3D is a rotor viscous code for three-dimensional 
turbomachinery flows which can only solve for flow through an isolated 
blade row. It solves the thin-layer Navier-Stokes equations with an 
explicit finite difference technique. The thin-layer assumption allowed 


for the streamwise viscous terms to be neglected. All cross-channel 
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viscous terms are retained. A Baldwin-Lowmax algebraic turbulence 
model was used. The equations were solved by using second-order finite- 
differencing spatial discretization. A multistage Runge-Kutta scheme was 
then applied to time march the solution. A more detailed description of 
RVC3D can be found in [Ref. 6: 1992]. 

RVC3D was used by Greco to simulate flow through the first-stage 
stator and rotor separately. The modified stator described above was run 
with identical flow conditions as specified by Greco [Ref. 4: 1995]. The 
output 1s presented in Appendix A. | 


a SWIFT 


The SWIFT program ts a three-dimensiona] thin-layer Navier- 
Stokes program for multi-blade row turbomachinery flows. It also used 
finite- difference formulations with an explicit multi-stage Runge-Kutta 
scheme that implemented variable time-step and implicit residual 
smoothing. The program had the capability to run multi-blocked grids 
containing C-grid around blades, H-grids upstream, O-grids in the hub and 
tip clearance regions, and mixing-planes between blade rows. The 
program contained the Baldwin-Lomax and Cebeci-Smith algebraic 
turbulence models and the Wilcox k-omega turbulence model. A full 
manual for the program had not been released at the time of the present 
study. For details on running the program see [Ref. 9: 1997]. 

The program used four input files. The fort.1 file contained the 
grid coordinates. The “fort.2” file was used for restarts and was a copy of 
the solution file produced by SWIFT. The “fort.10” file contained 
information on grid interaction and rotation. The namelist file 
“muliswft.in” contained the initial conditions for each grid, as well as 
program flags. The program produced two files. The “fort.3” file was the 
solution file and was copied to “fort.2” for restarts. The “swift.out” file 


was an ASCII file that contained output information on the run. 


2 


The SWIFT program was initialized with calculated values from the 
inlet pressure and Cobra probe exit-plane survey. Momentum averaging 
was selected to pass information between the grids. The input lists for the 
runs are included in Appendix B. The run was started with a Courant 
number (CFL) of 2.0 and a first-order artificial viscosity of 1.0, for the 
first 100 iterations. The program was then restarted with a first-order 
artificial viscosity of 0.5 for another 100 iterations. The first-order 
artificial viscosity was then further dropped to 0.25 for the next 100 
iterations. Zero was used for the first-order artificial viscosity for the 
next 100 iterations. The CFL was then change to 4.0 and the number of 
stages for the Runge-Kutta scheme was changed from two to four. The 


program was then run to 5000 iterations. 
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IV. RESULTS AND DISCUSSION 


A. COBRA PROBE MEASUREMENTS 


Four first-stage rotor exit Surveys were conducted using a three- 
hole (Cobra) pressure probe. Reference rotational speeds ranged from a 
low of 4863 rpm to a high of 4869 rpm. SSME and TTR data were 
collected for each survey and are presented in Appendix D and E. 
Appendix F contains the VEL_PRFL.VI data output as well as graphs of 
Mach number and swirl angle vs. radial position for each individual 
survey. Figures 23 and 24 show the combined graphs of Mach number and 
swirl angle vs. radial position for the four surveys. 

Both the Mach number and swirl angle graphs showed good 
consistency in the measured results and did not show any dependency on 
circumferential location. The complex structure exhibited by both graphs 
suggested the existence of secondary flow effects. The most plausible 
explanation of these structures are tip-leakage vortices generated by the 
rotor tip gap. This gap was measured, while rotating, to be .045 inches 
after resurfacing of the inner wall. This differs from the measured value 
of .020 inches, measured while rotating, reported in Southward [Ref. 5: 
1998]. The inner casing wall was resurfaced due to damage caused by 
bearing failure. 

The distribution of both plots were similar to that of Southward’s 
[Ref. 5: 1998] measurements. A more pronounced “s” shape was observed 
in the plot of Mach number vs. radial position compared to Southward’s 
plot. Minimum Mach number was observed at 89% (vs. 80%) span while 
maximum Mach number occurred at 65% (vs. 20%) span. This can be 
attributed to a larger vortex being generated by the larger tip clearance 


gap. The shape of the swirl angle plots were almost identical, however 


Zo 


the angles measured were approximately -5° lower throughout when 


compared to Southward. 
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Figure 24. Combined Cobra Probe Measurements of Swirl Angle 
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B. LASER DOPPLER VELOCIMETRY MEASUREMENTS 


The LDV surveys were conducted at two axial locations in the rotor 
tip-gap region. The majority of the surveys concentrated on the forward 
hole, located at —0.16c,. The other survey position was located at 0.35c,. 
The reference rotational speed ranged from 4781 rpm to 4862 rpm. Data 
were taken at three depths per axial position. The three depths, from the 
end-wall, were located at 0.019 inches (98% span), 0.069 inches (93% 
span) and 0.119 inches (88% span). Refer to Southward [Ref. 5: 1998] for 
a detailed description of the setup and operation of the LDV. 

Repeatability of the data sets was not good. The tangential flow 
velocity and angles varied widely. Only the forward position’s innermost 
depth had some repeatable data. This set was the axial component of flow 
velocity. Three of the four measurements varied only 4%. All other data 
sets followed no pattern. Figures 25 through 36 show the results of the 
plotted data. Appendix G contains the window-averaged data with 
parameters for each survey. With the exception of one forward inner LDV 
survey, the tangential laser beam was shifted by -1OMHz. For all runs the 
axial beam was not shifted. Both processor filter settings were set to 5- 
30MHz for all runs. 

Both frequency-shifters were found to be faulty and were sent out 
to be repaired. Upon return and testing of the repaired frequency-shifters, 
the PHASE optical setup was found to be in error. An input of —10 had 
been used as the setup for the tangential frequency shifting. This number 
should have been a positive 10. Before any other measurements could be 
preformed, the axial compressor drive motor developed a vibration and 
was sent out for repair. This precluded any further measurements on the 


TTR. 
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Figure 25. LDV Velocity Ratios for -0.16c,; and 88% Span 
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Figure 26. LDV Absolute Flow Angle for -0.16c, and 88% Span 
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Figure 27. LDV Velocity Ratios for -0.16c, and 93% Span 
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Figure 28. LDV Absolute Flow Angle for -0.16c, and 93% Span 
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Figure 29. LDV Velocity Ratios for -0.16c,; and 98% Span 
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Figure 30. LDV Absolute Flow Angle for —0.16c; and 98% Span 
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Figure 31. LDV Velocity Ratios for 0.35c, and 88% Span 
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Figure 32. LDV Absolute Flow Angle for 0.35c,; and 88% Span 
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Figure 33. LDV Velocity Ratios for 0.35c, and 93% Span 
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Figure 34. LDV Absolute Flow Angle for 0.35c, and 93% Span 
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Figure 35. LDV Velocity Ratios for 0.35c, and 98% Span 
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Figure 36. LDV Absolute Flow Angle for 0.35c, and 98% Span 
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c. NUMERICAL RESULTS 


i Modified Stator Model 


Once the grid was developed to model the step in the outer casing, 
both RVC3D and SWIFT were run to 5000 iterations and the data 
compared. The residuals for the RVC3D had not reached convergence but 
the SWIFT model converged after 2500 iterations. The residuals for the 
RVC3D run dropped over three orders of magnitude while the SWIFT had 
stabilized at 2.8 orders of magnitude. Plots of the residuals and 


coefficient of pressure at mid-span can be found in Appendix A. 


a. Combined Model 

The combined first-stage grid was run out to 5000 iterations. The 
residuals reduced by two orders of magnitude before converging at 3500 
iterations (Figure 37). The exit plane of the rotor showed very little grid 


dependency as can be seen in Figures 40 through 46. 
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Figure 37. Residual of Combined First-Stage Stator and Rotor 
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The SWIFT program produced a momentum-averaged output of the rotor 
exit plane. These data are included in Appendix H and were used to plot 


the Mach number and swirl angle vs. radial position (Figures 38 and 39). 
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Figure 38. CFD Rotor Exit Plane Averaged Mach No. 
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Figure 39. CFD Rotor Exit Plane Averaged Swirl Angle 
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Figure 40 is a composite plot of blade surface pressure distributions 
for the stage. The midspan Mach number distribution is plotted between 
the blades and the exit plane Mach number distribution is shown as 
contour lines. The high-pressure stagnation region is evident on the 
leading edge of the stator blade. Figure 41 shows the stage midspan Mach 
number distribution, which shows increasing magnitudes as the flow 
accelerates through the turbine. In these, figures the contour lines show a 
distinctly circular pattern of Mach number gradients downstream of the 
rotor tip. This is indicative of a vortex being shed by the rotor. In Figure 
42, tip leakage can be observed between the rotor tip and the end-wall. 
The flow in the step-region decreases in Mach number as it passes the 
step. This result would be expected. In Figures 43 through 45, the static 
pressure is shown on the rotor wall. In these figures, it can be seen that 
the pressure on the rotor tip drops as the flow speeds up and leaks 
between the tip and the end-wall. The tip leakage would support the 


generation of a trailing vortex. 
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Figure 40. Combined First-Stage (Static Pressure on Walls, Mach Number at Mid 
Section and Mach Contours at Exit Plane) 
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Figure 41. Combined First-Stage (Mach Number at Mid Section, k=24) 


Mach No. 


00 02 05 O07 1.0 





Figure 42. Combined First-Stage (Mach Number at Mid Section, k=50) 
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Figure 44. Combined First-Stage (Stagnation Pressure Contours at Exit Plane) 
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Figure 45. Combined First-Stage (Static Pressure Contours at Exit Plane) 
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D. EXPERIMENTAL VS. NUMERICAL COMPARISON 


Momentum-averaged computational results for the rotor exit plane 
were plotted for comparison with the Cobra probe measurements (Figures 
46 and 47). The predicted Mach number distribution in Figure 46 showed 
good agreement with the measured data, particularly in the tip region, 
where the radial location of the tip leakage vortex was correctly 
predicted. The computational results did not show the distinct “s” shape 
in the mid-section that the Cobra probe data exhibited, but the average 
magnitudes were fairly close. The grid line spacing in the mid-section of 
the rotor was large and this could be a possible cause for the departure. 

Figure 47 showed the swirl angle results. In the end-wall region, 
the plots overlapped. The numerical model did not follow the Cobra 
probe measurements from 60% span down to the hub. This departure can 
be attributed to poor modeling of the hub region. The step inward of the 
rotor hub, in the test turbine, was not modeled. This was the greatest 
source of disagreement between the experimental and computational 


results in the hub region. 
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Figure 46. CFD and Combined Cobra Probe Measurements of Mach No. 
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Figure 47. CFD and Combined Cobra Probe Measurements of Swirl Angle 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. COBRA PROBE MEASUREMENTS 


The Cobra probe measurements were successfully completed. The 
results showed that the inlet struts (of which there were 16) had little 
effect on the downstream rotor exit flow. The replacement of the soft 
plastic tubing on the probe increased the reliability and repeatability of 
the measurements. Redesign of the Cobra probe mount is needed in order 
to protect the probe tip while being installed and removed. A method for 
a more precise axial alignment of the probe is needed, to eliminate the 
inaccuracy introduced by only visually aligning the probe tip with the 


axial direction of the turbine. 


B. LDV MEASUREMENTS 


The LDV measurements were considered to be unsuccessful. The 
non-repeatability in the LDV measurements can be attributed to the noise 
that had developed in the frequency shifters. Replacement of these 
frequency shifters should eliminate this problem. The modification to the 
seeding wand did increase data rate acquisition between blade passages. 
The correction of the PHASE software optical setup should also lead to 
better correlation between frequency-shifted data and non-shifted data. 
Backscatter is still a problem. A field stop should be installed as 
Southward recommended [Ref. 5: 1998]. The refurbishment of the axial 
compressor drive motor should help stabilize the TTR rotational speed, 


leading to smaller variations in reference rpm between successive runs. 


c. COMPUTATIONAL FLUID DYNAMICS 


The initial computational fluid dynamics modeling results were 
successful. Good agreement between the Cobra probe measurements and 


the computational results was obtained. Insight into the flow around the 
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rotor tip was also gained. The stator and rotor grids both need 
refinement. Addition of more span-wise points (increase the dimension of 
k) is needed for both grids. The stator grid needs greater clustering of 
points in the axial direction, aft of the stator’s trailing edge. The rotor 
grid needs an increase in the density of points in the j dimension in order 
to better model the step in the end-wall casing. Investigation of the other 
turbulence models available in the SWIFT program could lead to better 


agreement between the experimental and computational results. 
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APPENDIX A. MODIFIED STATOR 


A. TCGRID INPUT FILE FOR MODIFIED STATOR 


&naml merid=0 im=301 jm=31 km=45 it1=95 icap=12 é&end 

&nam2 nle=19 nte=16 dsle=0.010 dste=0.0050 dshub=0.0001 
adstip=0.0001 dswte=-0.005 dswex=0.03  dsthr=120 deman=0-0061 
dsmax=0.02 rcorn=.098 &end 

&nam3 iterm=150 idbg=0 0000000 aabb=1.0 &end 

Enamaezbc=0,.0000 070000 272500 02000050. 0000 2227500 
rpc=4.-0788 4.0788 4.07638 5271430 5.1480 5.14680 fend 

‘new data style with z,th,r format SSME HPFTP ** COURSE GRID **' 

58 56 
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Drer=.9 cmutm=14..0 jedge=15 kedgqe=22 a11tin=2 abirn=0.005 
dbilt=0.005 &end 

Emil Gnixyd=0 oar=02 209251 102=151] nj0=1] nko=3° jJO=heko-13 17 21e8end 


C, SWIFT INPUT FOR MODIFIED STATOR 


'SSME HPTFP ADT FIRST STAGE STATOR' 

&ni2 nstg=4 c£il=4.0 aviscl=02 0 7avisc2=125 aviscc4=1. 50 oa 
eps=0.75 epi=1.0 ep) =1.0 7epk=1.0 pek=0-15 itmax=5 000 
Tvet=1 neeses=1 teims=.2 moec-0 cend 

SEnUse bein=3 abeex=-4 isyme=-O,8res=1 reser © 
Presto= il Bwoin—0 cend 

fni4 aqeom=1 ga=).4 om=0 .0 prae=07665 eroee0r 0 &end 

&EnlS5 alte=2 renr—-e624075.8 prar=0 27 tw=lmOmem cowr=.66666- 
Dme@r=.9 cmuim=14.0 jedge=15 kedgh=22 kedgt=22 ialein—2 
abih=0.005 "dblt=0.00S frough=4.0 mreur—5 send 


&nl6 oar=0.0 mioe=3 igqav=1 é&end 


row PO Mx Mt Mr TO 
0 170000 O22 126 oO 0:0'6 0 0.0000 1. 0G 
al 0.9850 m2. > 0 -0.6720 0.0000 1. 00CG0 


D. MODIFIED STATOR FORT.10 FILE FOR SWIFT 


grid type im jm km al ia aa nmin nex nhub ntip noalr row om omh omt 
a Z 135 oon 7 dz 62 0 Sous 99 0 0 0 a ie 2 O'. 
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EK. MODIFIED STATOR RESULTS 
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Figure A.1. Modified Stator Grid (301x31x45) 


ay 


-20 





-25 
0 0.2 0.4 0.6 0.8 1 


Figure A.2. RVC3D Solution of Coefficient of Pressure at Mid-Span for Modified Stator 





Figure A.3. SWIFT Solution of Coefficient of Pressure at Mid-Span for Modified Stator 
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Figure A.4. RVC3D Residuals for Modified Stator 
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Figure A.5. SWIFT Residuals for Modified Stator 
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APPENDIX B. MGRID PROGRAM 


The MGRID program is a Fortran program for repeating grid files 


circumferentially. The output can then be viewed in NASA’s FAST 


program. 

Cc — om em em em oe me me mm ma wm ww ww we we we sw oe ee ee 8 
c Prog. reads an unformatted multigrid fort.1 file and produce 

C a multigrid file of n blades. Change n to the desired number of 
e blades to be shown and change h to the number of blades on the 

C turbine or stator. 

Fc enc es pg a ee Oe 


tb St Se SE SE 


te 4b 4 SE SE 


real (300,100,100) vy (3007 h00 200) 4 zi S00 3100 200) 
realryyt200, 100;5100)).2 21300, 100, 100) 
integer ni(200),nj(200),nk(200),h,1,n 
open (unit=21, file='mgrid.dat',status='unknown' ) 
n=3 
f= 50 
read(1) ngrid 
Deimta * grid 
Read(1) (nieGigrid) 1) (igeid). nk (iorid) .1Orid—1nourd) 
DEime ten (igrad), ny (igrid) maki iorid) 1646. 1d=),ngric) 
=n *norid 
write(521)1 
Drie. 9 ok 
weree (S21) ((nitioqmidy, ni (torrdj ukiigrid), m=), n).1orid— maa) 
Disime.* > ((maigrid)- ni (iorid) zok(tgqrid) m—-1,0)4 19grild-l) nerd) 
dont grid=1) norid 
mead (1h) 0 (alah 1c), 


Paden igerd) yo j=longjtigrid) 4 keteank( tora). 
VACA aac) 
P=lerd (19rd) ) jal, ay torrd)) k=) ak(i¢gelain: 
CUGZ ae je, 
=, Win Loewe) ) 4 jay (a Or) ye ie Peete cpr eciyey 
pid2=acos(0.0) 
do m=1,n 
bang=(m-2)'*4..*pid2/h—2 2*pid27h 
do K=1, nk (grid) 
Cot = lint chia) 
d6-j- Ly torid) 
Vy (1,3-,K)=y¥(2, 32k) *cos(bang)+z(i,7,k) “sinibang) 
zZz{l,),K)=-V (isd, kK) *Sin (bang) 4+z2(@, 1, kK) * eos (eanc. 
enddo 
enddo 
enddo 
Write (S24 Gi (agg, ele, 
i=1,ni(igrid) )7j=1-nj (i90rid)) k=l nko. 
GOV (aaa 
i=], ni (worid)),j=1;nj (igrid)) k=) neaicria 
CGlzz ia tee 
t=1,ni (igrid) )74=1 nj (29rid) je k=l nk 
enddo 
enddo 
stop 
end 
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APPENDIX C. COMBINED NUMERICAL INPUT FILES 


A. NEW TCGRID INPUT FOR COMBINED STATOR AND ROTOR 


a Stator.in File for New TCGRID 


&naml merid=0 ‘im=135 jm=31 km=57 itl=12 icap=12 igcelt=0 
Kmt=13 jJmt=135 ég&end 

&nam2 nle=19 nte=16 dasle=0.010 dste=0.0050 dshub=0.0001 
aistip=0.,000l=dswte=-0 005 7 aswex=0.03 dasthr=l)s.07dsmin=0_001 
asmax=0,005 reorn=.0938 cltip=.045 €Gdscit=.0002) tend 

RMeams. Lteerm=150 a2G@bgq=0 0°0CO 0 °070702 050 aabb=1 .0 send 

Srem4. 2zpc=0.0000 0.0000 T5000 7020000 0.0000 175000 
roc=4.0788 42,0788 4.0788 SS. les07S 14860 527024725 “send 

&namS iswift=1 dslap=.005 é&end 

“new data style with 2,2h,r format SSME HPFTP ** COURSE GRID *** 

Se Boo 


-0.9334116 -0.8467470 -0-7600824 =0 2673534079 =O. 5567434 
=O. 0.000 84 -0.4134042 -VzSs25 7396 =O 2400750 =0 21534005 
-6.673590le-02 1-9928699e-02° 86,1 066032 Use 3 2676 Osz 799324 
Ui 6 66015 9 O-4552 715 O25 399561 Oy 6:26 60107 0.6927030 
C2 Os 79 OF 7360452 027577163 0. 7837138 Vas 027210 
Oi S352 1 O58 5350533 G2879 0605 0.2007 Sid Oo 267292 
0.9527364 Oe Se o-4 37 1.004741 1: ,025.055 1 061082 
i208 27 43 eed 51-4 8 1.216443 La 239866 12655918 
1287020 eo 0 7 3.0 Ue sraia.o, 0 l.2o654 9 Pe See 228 
Oro 1 37 1743.00.51 aod oO De Odeiaio2 1.620448 
1622669 1.624487 P2625 05'5 ive aoe 3 met 6 2 24-1 
(1.629662 ewe Sao bd ifs 0000 

Seo 294 Seve 32 2550925 32.0640 6 Bu-68 6 206.01 
32902798 Bio 227.46 Be soo Jaga 6246 Sano eeee.e) i 
32986304 52999689 4.012004 4.023261 4.033448 
42042595 4.050704 Oe 115 4063811 40067732 
4.068662 A069 97 9 a2 T0999 4.072147 472057 5206 
4.074018 4.074761 a7 075573 4.076187 42/0:76830 
4.077384 42077650 A107 2 16 4.0771 89°33 42078711 
4.078780 4.078800 4.078800 4.078800 42107 6:5.00 
4.078800 4.078800 4.078800 4.078800 4.078800 
4.078800 4.078800 4.078800 4.078800 4.078800 
420758600 4.078200 4.078800 4.078800 4.078880 
4.078800 4.078800 4.078800 

O19 35.5 4416 -0.8467470 -0.7600824 =0..673407 9 =-0-5867434 
=O 5000/87 -0.4134042 -0.3267396 =O 24.0075.0 -~QO2 2544005 
“oo cou lLe-OZ 1.9928699e-02 0.1066032 OP os 2678 062799324 
023665069 04562 71.5 Oo 2993 61 076266007 Oms6:2 75.5 5.0 
O2711473'4 Ji so5 819 Ue 592607 Ol Grs.ci 9-9 Ose0G e797 
Gos 05002 0.8542116 036777439 Oe 2011673 0.9244521 
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PPP PP oO 


wun nnn ongo sol ll ll 


o) @) (e) fe) (©) fee me 


rFPrPrrRFooosoeo0oee#eoeo0nouoed0e#9W0doe9°coao 902 OF FF F FF 


ooo 0) 6 
.085674 
.287020 
~406-13 7 
O22 oom 
56:2 9'662 


.148000 
mag O00 
.148000 
.148000 
wie 713 
we. 2-0.155 
7 lyse Ast yas 
moos. 2 
.034636 
.024538 
.026436 
oU9G3 07 


oS 52 
m25353 8 
meei2453 
poo 042 
70 563.3 0 
7D 618 
oo 0 90'S 
P14 8192 
.645480 
-542768 
£69357 
-448272 
pad ZI23 
754.0108 
so043.95 
nt O7 107 
moO l 820 
aoe oo 3 
2025245 
mile? 958 
leo 19 
meow 104 


-091840 
-083090 
.064230 
aos 9170 
mera) 30 
2011340 
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201027 30 
.007360 


PrPerPrPrP Oo 
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OS 637 
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204236 
wo 500 
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-440055 
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E22 ee 
27624937 
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we2e 5 21 


wog2470 
wo 9820 
.058240 
Co 3 40 
mo 20 120 
moo.9 9.0: 
wo01ls60 
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-006800 


PPP PP RP 


nono uno nun noun wo WM 


RPrerPrRreooooeoeqecjceooooo0o090o 02 0 FF FF F 


~OO05068 
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190000 
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sta 8000 
48:00:10 
la Do. 6 
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ei 8] he ee ag A 
047565 
wO/Ze-87 3 
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Oy ls 
2 09'S 87.5 


a2o7 104 
~196019 
doo 
205245 
woe oS 3 
~509820 
oO 107 
e043 95 
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n0:9 500 
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Dre Combined Stator, Rotor and Tip Input File for SWIFT 
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TTR Data (2 of 2 


[Date [File Name Water inlet [Water Outlet swirl angie [throttle pos atm press (in Hg 
[227706 537.279247 1546. a 


2/27/88 567.695292 : 541.283492 
12/27/98 568.06266 


W398 559.314181 
WV98 §65.425434 
13/5/98 


535.2998 
543.453241 
$42.765572 
575.227595 |565.74395 
0310fwd1.10d]566.842572 1545.157558 
031Ofwd1.20d }568.740734 


579.871173 
574.826145 
577.756429 


13.3490713 
12.9035284 
13.1626058 
13.3653672 
13.470313 
11.5888709 
19. 143854 
13.1111743 
13.0000727 
11.0596747 
12.6445781 
12.6560713 
11.6723909 
13.1911188 
13. 1288255 
13.1379684 
13.0507524 
10.9990266 
11,.4039859 
13.186107 
13.2393709 
13.150341 
13. 1052403 
13. 183305 
13.1090442 
10.6285508 
12.4414992 
12.409931 
12.4493037 
12.36662 


281.47113 


10.6347656 
10.4589844 
10.1953125 


13.6230469 
14.3261719 





Table E.2. TTR Data (2 of 2) 
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APPENDIX F. COBRA PROBE MEASUREMENTS 


Radial Position vs. Rotor Exit Mach No. Position 1 (05 Mar 98 run2) 
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Figure F.2. Cobra Probe Mach Number: Position 2 (27 Feb 98) 
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Radia! Position vs. Rotor Exit Mach No. Position 2 (05 Mar 98 runt) 
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Radia! Position vs. Rotor Exit Mach No. Positon 3 (3 Mar 98) 
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Figure F.4. Cobra Probe Mach Number: Position 3 (03 Mar 98) 
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Radial Position vs. Rotor Exit Swirl Angle Position 1 (05 Mar 98 run2) 
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Figure F.5. Cobra Probe Swirl Angle: Position 1 (05 Mar 98) 


Radial Position vs. Rotor Exit Swirl Angle Position 2 (27 Feb 98) 
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Figure F.6. Cobra Probe Swirl Angle: Position 2 (27 Feb 98) 
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Radial Position vs. Rotor Exit Swirl Angle Position 2 (05 Mar 98 run1) 
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Radial Position vs. Rotor Exit Swirl Angle Position 3 (3 Mar 98) 
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Figure F.8. Cobra Probe Swirl Angle: Position 3 (03 Mar 98) 
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SSME HPFTP = ——— — Exit Velocity Survey (Position: 1 } 


03/05/9Brun2 


-13.5083333 | 135.324333 117.622 14.314 4.87866667 1.02234375 | 5083.10578 | 0.02179101 | 0.10800669 | 0.24293141 | 0.77488278 








































































































































































































































































































































































































































-13.597 135 .334 117.577333 | 15.6843333 2.186 : 5095.76702 | 0.03106986 } 0.13004374 | 0.29327708 § 0.76308291 
-13.7456667] 135.298667 J 117.499333 | -1.45866667 0.94421875 | 5094.4584 } 0.01660851 | 0.09569708 } 0.21497178 § 0.75357468 
-13.9376667] 135.296667 | 117.139667 | -11.4263333 0.90027344 | 5093.65466 f 0.00911742 | 0.07950013 | 0.17833214 | 0.74356358 
-14.0806667§ 135.35 7667 117.486 | -10.7393333 25.5761719 | 0.86121094 | 5098.805 [ 0.01043872 | 0.08212556 | 0.18426078 | 0.74328613 
°14.3176667] 135.292667 | 115.561333 -4.6 36.4746094 | 0.81238281 | 5067.18342 | 0.01473653 | 0.09140133 | 0.20523868 | 0.75298255 
-14.4806667 ] 135.305667 | 115.445667 | 4.17366667 43.68 16406 | 0.77820313 | 5065.72158 f 0.02366627 | 0.11251958 | 0.25320944 | 0.75750656 
-14.6786667 | 135.349667 | 115.453667 | 13.6586667 46.6699219 | 0.73425781 | 5080.03562 | 0.03434623 | 0.13746049 | 0.31031675 | 0.7575318 

135.283 115.187 | 19.3106667 47.9003906 | 0.69519531 | 5071.23684 | 0.04007645 | 0.14969492 | 0.33854263 | 0.75831455 
135 292333 | 114.891333 21.448 47.9003906 | 0.65125 | 5078.25542 | 0.04203499 | 0.15363627 | 0.34766885 ] 0.75923444 
135.310333 115.129 | 20.7946667 46.4941406 | 0.61707031 | 5086.83954 § 0.04210972 | 0.15378413 | 0.34801154 | 0.75770334 
135.289667 | 114.617333 | 18.3856667 45.3515625 | 0.57800781 | 5072.48702 { 0.03960361 | 0.14872443 | 0.33629802 | 0.75848645 
135.280333 | 113.886667 15.129 44.296875 | 0.53894531 | 5066.19702 | 0.03578388 | 0.14062369 | 0.31760009 | 0.7602812 
135.285333 ] 114.259667 12917 45.1757813 | 0.49988281 | 5068.65462 | 0.03310189 | 0.13467606 | 0.30391357 | 0.76036868 
135.347 114.174 11.65 49.1308594 | 0.46082031 | 5081.59216 | 0.03096771 § 0.12980828 | 0.29273694 | 0.76179933 
-16.1276667] 135.282667 | 114.393667 } 11.1243333 54.0527344 } 0.42175781 | 5086.48642 | 0.0312095 | 0.13036524 | 0.29401465 | 0.76021774 

-16.322 135.273 9 113.890667 | 11.1186667 59.4140625 } 0.38269531 } 5074.93686 | 0.03233472 | 0.13293909 | 0.2999229 | 0.75915324 
-16.4993333] 135.278333 | 113.706667 | 11.4046667 63.1054688 } 0.333867 19 | 5075.44406 | 0.03379102 | 0.13622326 | 0.30747067 | 0.75758486 

-16.666 135.34 113.572667 } 11.6173333 66.0058594 | 0.28992188 | 5081.7052 | 0.03476117 | 0.13837952 | 0.31243184 | 0.75660237 

-16.805 135.283 113.494 11.58 67.2363281 | 0.25574219 | 5085.68356 | 0.03564487 | 0.14032044 | 0.31690143 | 0.75523497 

-16.943 135 285 113.584667 | 11.5706667 67.2363281 | 0.20691406 | 5082.95952 | 0.03642673 | 0.14201859 | 0.320815 | 0.75384479 
-17.1323333 | 135.266333 113,092 11.889 66.5332031 | 0.16785156 | 507138678 | 0.03695992 | 0.14316609 | 0.32346118 | 0.75427246 
-17.2653333] 135.273333 112.816 | 12.8043333 64.7753906 | 0.12878906 | 5076.95566 § 0.03792013 | 0.14521049 | 0.32817896 | 0.75470358 


Table F.1. Cobra Probe Data: Position 1 (05 Mar 98) 


SSME HPFTP ATD First-Stage Rotor Exit Velocity Survey (Position: 2 } 


a a a 
2/27/98 


























































































































































































135.902667 | 131.213667 | 26.6313333 } 15.6113333 -3.1640625 § 0.96375 | 5091.86712 | 0.02522352 | 0.11626043 | 0.26174116 | 0.76932699 
135.900667 ] 130.427 12.098 4.81166667 -0.96679688] 0.92957031 | 5073.1986 ] 0.01723721 | 0.09716387 | 0.21829792 } 0.75577988 
135.870667 | 130.059 {| 1.64066667 | -3.32766667 10.0195313 | 0.885625 | 5081.24696 | 0.01203315 | 0.08545021 J 0.1917739 | 0.74344587 













































































































































































































































































































































135.857667 | 129.759667 | 0.90566667 -4.361 18.28125 | 0.8465625 | 5095.89802 J 0.01276773 } 0.08703199 | 0.1953507 } 0.74196725 

135.858667 | 129.637667 7.631 0.929 27.8613281 | 0.80261719 | 5084.18192 | 0.01597481 } 0.09422996 | 0.21164632 | 0.75089527 

: 135.816333 } 128.774667 | 17.2643333 6.685 37.0019531 | 0.75867188 | 5099.86464 | 0.02459455 | 0.11475119 | 0.25829772 | 0.75813111 
-3.86966667] 135.81 7667 | 128.539667 | 26.7036667 12.048 42.8027344 | 0.71472656 | 5099.27474 | 0.0333112 1 0.13514733 | 0.30499681 | 0.76130391 
-4.23966667 | 135.817667 | 127.808333 30.429 13.896 449121094 | 0.67566406 | 5087.3482 | 0.03723167 | 0.14374757 | 0.32480261 | 0.76234826 
-4.802 135.815 127.396 29.605 13.366 46.494 1406 | 0.62683594 | 5099.43264 | 0.03659115 | 0.14237337 | 0.321633 | 0.76275887 
-§.19233333] 135.786333 127.1 28.73 12.8416667 46.8457031 | 0.61707031 | 5094.9267 | 0.03584014 | 0.14074626 | 0.31788251 | 0.76304759 
-§.51033333] 135.776333 | 126.825333 | 26.2176667 j 11.2693333 46.3183594 | 0.573125 | 5096.00538 | 0.03388747 | 0.13643881 | 0.3079664 | 0.76245155 
-5.85566667] 135.813667 | 126.968 24.784 10.3706667 44.5605469 | 0.52917969 | 5101.3642 | 0.03275546 | 0.13389356 } 0.30211543 § 0.76175221 
-6.342666671 135.762667 | 125.945667 | 24.0493333 | 9.48966667 41.6601563 | 0.49011719 | 5096.79814 | 0.03310664 | 0.13468678 | 0.30393822 | 0.76149741 
-6.86866667 | 135.773667 | 126.248333 23.629 9.085 45.2636719 | 0.4559375 | 5103.04602 | 0.03306308 | 0.13458855 | 0.30371245 | 0.76058864 
-7.13166667 | 135.756667 | 125.033333 | 22.5143333 | 8.00133333 48.8671875 | 0.41199219 | 5101.55558 | 0.03305661 | 0.13457395 | 0.30367891 ] 0.76046116 
-7.824 135.74 124.829333 | 20.8266667 6.912 54.4921875 | 0.37292969 | 5095.50692 | 0.03176578 | 0.13164148 J 0.29694347 { 0.76016513 
135.726333 124.64 19.7076667 | 6.19166667 59.3261719 | 0.32898438 j 5106.47568 | 0.03090694 |} 0.12966809 | 0.2924154 § 0.76001421 

135.726667 | 124.202 | 19.8946667 § 5.75233333 63.3691406 | 0.28992188 | 5104.37434 | 0.03229971 | 0.13285947 | 0.29974004 | 0.7588026 

135.707 | 123.716333 } 20.6176667 } 5.54933333 66.6210938 | 0.25085938 | 5105.73712 } 0.03432982 | 0.13742403 | 0.31023286 | 0.75751609 

135.723667 | 123.378667 20.918 §.24733333 68.0273438 j 0.20691406 | 5109.26694 § 0.03564257 | 0.14031544 | 0.3168899 | 0.756477 

135.693333 | 122.864667 | 20.7806667 5.035 67.5 0.17273438 | 5094.42762 | 0.03580369 | 0.14066686 | 0.31769955 | 0.75677385 

135.727 122.505333 20.545 4.807 65.2148438 | 0.12390625 | 5110.6022 | 0.03577178 | 0.14059731 | 0.31753931 | 0.75705712 

135.695333 | 120.594 19.0943333 | 3.19166667 64.86328 13 | 0.09460938 § 5097.76176 | 0.03611768 | 0.14134953 | 0.31927274 | 0.75723903 

135.623333 | 119.819667 | 15.5563333 5.18554688 | 0.98816406 | 5110.87922 | 0.02090114 | 0.10586861 | 0.23806731 | 0.77627385 

135.640333 ¥ 119.267667 -13.2714844] 1.03210938 | 5099.18338 | -0.0048619 | 0.06431108 | 0.14410226 | 0.74212922 

135.650333 | 119.446667 15.1171875 | 0.86121094 | 5100.47522 | 0.01210678 | 0.08560742 | 0.19212933 | 0.7406485 

-13.2856667] 135.604667 118.761 -12.8436667 23.7304688 : 0.01374271 | 0.08917389 | 0.20019644 | 0.74633206 
-13.469 135.652 | 119.093667 -6.90166667 32.60742 19 | 0.78308594 | 5105.9087 | 0.01847321 § 0.10007506 | 0.22490368 | 0.75489249 
-13.628 135.608 ] 118.467333 -0.35666667 40.3417969 | 0.73425781 | 5097.35864 } 0.02861345 J 0.12432118 | 0.28016413 | 0.76091613 


Table F.2. Cobra Probe Data: Position 2 (27 Feb 98) 
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a — ATD First-Staqge Rotor Exit Velocity Survey (Position: 2 ) 


fs ea i fee | Nedaees cee || 
sacs 


125.732333 | 24.6833333 : 5.9765625 § 1.01257813§ 5059.3154 | 0.02156441 | 0.10746184 | 0.24169155 | 0.77902069 































































































































































































































































































































































































































135.568 125.97 1667 1.49414063 | 0.98328 125 | 5091.66108 | 0.02396822 | 0.11324589 | 0.25486506 | 0.76397855 

135.546333 123.227 5.36132813 § 0.93933594 | 5053.43842 | 0.01503659 | 0.09208169 | 0.20677943 | 0.75350376 

135.539 122.684 12.1289063 | 0.90515625 | 5084.4813 | 0.01130248] O. 0.18828595 | 0.74222336 

-9.14066667] 135.531667 | 121.344667 | -6. 21.5332031 | 0.86609375 | 5074.92748 | 0.0110441 | 0.08336915 I 0.18707032 | 0.74264109 
-9.65733333] 135.523333 | 121.504333 | - 30.0585938 | 0.82214844 | 5081.5961 | 0.01335914 | 0.08832577 | 0.19827737 } 0.74841532 
-10.6623333] 135.527333 | 119.668667 } 4. 38.8476563 | 0.78796875 | 5072.62446 | 0.01984704 | 0.10334406 | 0.2323283 | 0.75668076 
-11.7093333] 135.506333 | 119.165933 | 12.9423333 f 6. 44.296875 | 0.74890625 | 5082.0204 | 0.02933776 | 0.12602098 J 0.28405609 | 0.75812331 
“12.401 135.506 | 117.970333 19.247 3.05666667 46.3183594 | 0.70984375 | 5067.04526 | 0.03707005 | 0.14340202 I 0.32400543 | 0.75842214 
“13.171 135.49 117.319 22.351 4.413 47.3730469 | 0.68054688 | 5093.2534 | 0.04071048 | 0.15098474 } 0.34152738 | 0.75889505 
-13.664 135.498 | 116.884333 | 22.2273333 | 4.01733333 46.8457031 | 0.63171875 | 5073.9534 | 0.04129317 | 0.15215842 § 0.34424492 | 0.75848111 
-14.237 135.467 116.401333 19.871 2.57433333 46.40625 | 0.59753906 | 5086.03586 | 0.038938134 | 0.14826572 | 0.33523742 | 0.75845646 
135.531667 | 116.274333 | 17.1693333 45.8789063 | 0.56335938 | 5086.25338 | 0.03635919 | 0.14187261 | 0.32047846 | 0.75912144 

135.462333 | 115.564333 14.476 45.7910156 | 0.52429688 | 5069.24266 | 0.03359981 | 0.13579526 | 0.30648643 | 0.76019397 

135.484333 | 115.354333 | 13.1523333 46.5820313 | 0.48523438 | 5086.6023 | 0.03319657 | 0.13488938 | 0.30440387 | 0.75898345 

135.472 | 114.630667 12.041 49.921875 | 0.45105469 | 5081.54368 | 0.0321 4007 | 0.13249603 | 0.29890542 | 0.75987183 

135.455333 | 114.450667 } 10.7516667 54.3164063 | 0.40710938 | 5096.72056 | 0.03102816 | 0.12994765 f 0.29305665 | 0.75977188 

135.430333 | 113.481667 § 9.91533333 54.140625 | 0.41199219 | 5072.80316 } 0.03126543 | 0.13049388 | 0.29430981 | 0.75954793 

135.432333 | 113.235 9.626 60.8203125 | 0.36804688 | 5078.61438 | 0.03101664 | 0.12992109 | 0.29299573 | 0.75985944 

135.441333 | 112.944333 | 10.3346667 64.8632813 | 0.32410156 | 5078.58372 | 0.08339466 | 0.13533492 | 0.30542805 | 0.75772753 

-17,6056667{ 135.431667 7 112.345 | 10.6193333 67.0605469 | 0.27527344 | 5070.93662 | 0.034747 46 | 0.13834924 | 0.31236213 | 0.75690773 
-17.75166674 135.432667 J 111.553333 | 10.4583333 68.1152344 | 0.23621094 | 5077.04222 { 0.03533408 0.31533501 | 0.75682587 
-17.9466667] 135.438667] 112.346 | 10.5706667 68.2910156 | 0.20203125 | 5069.56464 | 0.03572261 0.31729229 J 0.7553234 
-18.1003333] 135.417333 | 111.924333 67.1484375 | 0.16296875 | 5075.8252 | 0.03660255 | 0.14239794 | 0.32168964 | 0.75505934 
-18.2856667 | 135.458667 } 112.104667 65.9179688 | 0.12878906 | 5068.99806 | 0.0368061 | 0.14283593 | 0.32269968 | 0.75568503 
-18.5786667] 135.395667 | 111.608333 22.8515625 | 0.85144531 | 5068.11198 | 0.01177937 | 0.08491069 | 0.19055424 | 0.74386881 
-18.7416667] 135.427667 | 111.252333 13.0078125 | 0.90027344 | 5075.49472 | 0.01061007} 0.0824751 | 0.18505037 | 0.74346518 
-18.8593333] 135.394333 | 111.561333 6.24023438 | 0.92957031 | 5066.53114 | 0.01409507 | 0.08995884 | 0.20197297 | 0.75281418 


Table F.3. Cobra Probe Data: Position 2 (05 Mar 98) 


SSME HPFTP ATD First-Stage Rotor Exit Velocity Survey (Position: 3) 


eee ae eee a 
3/3/98 












































































































































































































































































































































































































-0.728 135.955 | 130.807333 12. Sane it -4.13085938 } 0.97351563] 5084.32692 0.25412278 | 0.76487569 
~1.666 135.976 | 129.416667 -0.61523438 | 0.93445313] 5084.37554 : 0.21187173] 0.75001395 
-2.43666667 | 135.955667 | 129.372333 “5. cece 8.701171875] 0.90027344 0.08810129} 0.19776948 
-3.07133333 | 135.934333 | 128.349333 -5.99166667 17.31445313] 0.86121094 0.09314106] 0.20917906 
*3.74233333 | 135.963333 | 127.444667 -1.51033333 28.56445313] 0.81726563 0.1001019 | 0.2249646 
-4.57966667 | 135.935667 | 127.075333 4.366666667 39.0234375 | 0.7684375 0.12017837 | 0.27068888 
-5.324 135.952 125.952 8.122 43.59375 |0.73914063 : 0.13989981 ] 0.31593247 
~6.11166667 | 135.937667 125.51 9.923666667 45.79101563] 0.69519531] 5101.54174 0.15027126] 0.33987612 
-6.54266667 | 135.933667] 125.148 9.922666667 47.109375 | 0.65613281] 5094.25062 0.14942118] 0.33790941 
-7.49933333 | 135.955333 | 123.188667 8.391333333 0.61707031 } 5085.97334 0.14522762 4 0.32821852 
-8.03866667 | 135.916667 | 122.513333 7.096 0.58777344 } 5082.58882 0.3192055 
-8.48633333 | 135.935333 | 121.964667 6.289666667 0.54871094] 5076.5553 0.31471092} 0.75981417 
-8.886 135.98 121.862 0.50964844] 5080.79188 0.31449794] 0.75952001 
-9.34866667 } 135.912667 121.08 0.47058594]} 5068.11598 0.31879651] 0.7582209 
-9.77133333 | 135.921333] 120.728 0.44128906 | 5081.59294 0.31570934 | 0.75826236 
-10.124 135.999 | 120.760667 0.40222656] 5078.9642 0.31614735 | 0.75619381 
-10.415 135.929 | 120.284667 58.8867 1875] 0.363 16406 | 5076.25462 0.31 1803 
-10.7413333 ]| 135.913333 | 119.632667 63.72070313] 0.324 10156] 5076.65522 0.31077603 
-10.9806667 | 135.918667 | 119.642333 66.88476563] 0.28503906 } 5078.62588 0.1409374 §0.31832294 
-11.2526667 | 135.929667 | 119.692333 68.81835938 5069.898638 0.14449094] 0.32651799 
-11.5173333 | 135.907333 | 119.101667 68.46679688 5074.2438 0.32787111 
“11.831 135.916 | 119.223667 68 .203125 5085.74354 0.32743156 
°12.168 135.918 | 118.514333 65.83007813 5077.51688 0.32892021 
-12.5116667 f 135.906667 | 118.267667 64.42382813 5073.32418 0.33406727 
-12.8626667 | 135.922667] 117.818 32.78320313 ; 5072.07914 ] 0. ; 0.23914475 
-43.1556667 | 135.897667] 117.654 20.0390625 | 0.85632813] 5085.60934 7 0.0156623 | 0.093511041 0.21001729) 0.7426827 
-13.433 135.859 116.959 -17.0006667 10.37109375} 0.89539063] 5080.71352 | 0.01398748] 0.08971858} 0.20142917] 0.74057527 


Table F.4. Cobra Probe Data: Position 3 (03 Mar 98) 


of 


APPENDIX G. LDV MEASUREMENT DATA 


Window Ave: _On 

Filter Setting: 5-3OMHZ Sint freqshifting 

V-turb__fAlpna___—*([90-Alpna_{Cuv_ | X_tneta— X_axial_| 
0 79.93) 88.69] 24.47 8.911 42.0 


47.974 





Table G.1. LDV Data: Forward Position, Inner Depth (10 Mar 98) 


9 












Ne ps= Onn 
Vu 30-Alpha 
1 56 L cE 









Table G.2. LDV Data: Forward Position, Middle Depth (10 Mar 98) 


94 


First-Stage Rotor LDV Data 
Date: 03/10/9 Axial Pos: -0.16ct 

N Ref: 4839 Span Pos: 98 

Vt 


8 
Filter Setungs:_5-G0MHz 
lU-mean_|v-mean [Uturb |v-turb [Alpha __]90-Alpha 
0 99.86 74.64 15.92 20.98 53.226 36.774 





Table G.3. LDV Data: Forward Position, Outer Depth (10 Mar 98) 


7) 












Filter Settings: 5-30MHz freq shitting: -10, 0 





Table G.4. LDV Data: Forward Position, Inner Depth (06 Apr 98) 
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: 16.17 
8.83 16.13 
8.21 14.19 
13.52 





Table G.S5. LDV Data: Forward Position, Middle Depth (06 Apr 98) 
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req shifting. -10.0 
[Usmean |v-mean [Uturo[v-turb Alpha [90-Alpha [Ow Ja 
0 140.87 74.8 8.3 5.47 : .96 





Table G.6. LDV Data: Forward Position, Outer Depth (06 Apr 98) 


98 













i. avn en eo ee 
V-mean 
6.5 126.9 82.04 






U-turb 
29 


V-turb = JAlpha 90-Alpha [Cuv —_—s EX _ theta 
9 14.08 8 


Table G.7. LDV Data: Center Position, Middle Depth (06 Apr 98) 
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Finer Sotungs: 6-0MHz____ ]reqshiftng: -10, 0 





Table G.8. LDV Data: Forward Position, Inner Depth (10 Apr 98) 


100 


First-Stage Rotor LDV Data Window Ave: on 


Axial Pos: -.16ct 
ret: Span Pos: 93% 
FiterSetings 5-SOMHz _________ireqshifting. -10. 0 
“tur 
0 141.43 73.54 14.61 14.05 





Table G.9. LDV Data: Forward Position, Middle Depth (10 Apr 98) 
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First-Stage Rotor LDV Data 


Span Pos: 98% 
Filter Settings: 5-30MHz 


— 





0 : 76.55 14.6 2.45 61.12 28.88 


0.092717 
0.092 
0.089786 


Table G.10. LDV Data: Forward Position, Outer Depth (10 Apr 98) 
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Date. 04/10/98 
Span Pos: 88% 





Table G.11. LDV Data: Center Position, Inner Depth (10 Apr 98) 


Theta [U-mean_[V-mean_ [U-turb_ [V-turb_ Alpha [90-Alpha_[Cuv_ | x_theta]x_axial_] 
7.8] 12358] 70.29 10.4 2.57) 60.367| 29.633 





Table G.12. LDV Data: Center Position, Middle Depth (10 Apr 98) 


First-Stage Rotor LDV Data 

020M ae ee 

V-mean Alpha 
f2 : 33 o 2 58.628 31.372 ; 





Table G.13. LDV Data: Center Position, Outer Depth (10 Apr 98) 
103 


= 020 in 


Filter Setting: 5-30MHz freq shifting: -10,0 
Theta___[U-mean [V-mean [U-turb_ |v-turb [Alpha __490-Aipna fou‘. 
0 : 77.58 6.13 19.38 59.534 6 : ‘ 





Table G.14. LDV Data: Forward Position, Inner Depth (01 Apr 98) 
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First Stage Rotor LDV Data Window Ave: On 
Date:_05/01/98 Axial Pos: -.16ct 











































































N_ Ref: 4808 Span Pos: 93% 


Vi= 794 Eoin S| 
Filter Setting: 5-30MHz frea shifting: -10. 0 





33.329 
29.805 
30.458 
29.677 
30.132 

30.69 
30.636 
30.764 
30.752 
31.006 
31.085 
31.058 
30.813 
30.935 
Bier 
31.262 
31.305 
31.303 
31.054 
31.218 
31.147 
31.037 
31.111 
31.164 
31.307 
30.909 

31.35 
31.215 
31-295 
30.852 
30.863 
30.646 
30.387 
30.422 
29.676 
30.454 
30.277 
29.986 
30.485 
29.673 
29,559 
29.722 

29.38 
28.678 
29.474 
28.768 
28.333 

28.82 
28.599 

28.27 
29.056 
29.075 
28.785 
29.276 
28.566 
29.046 
29.074 
29.184 
29.401 
29.532 
29.175 
29.179 
29.407 

29.21 
29.324 

29.39 
29.783 
29.547 


Table G.15. LDV Data: Forward Position, Middle Depth (O01 May 98) 
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Fitter Setting: 5-30MHz freq shifting: -10, 0 


“tur 
0 fi 87.52 49 





Table G.16. LDV Data: Forward Position, Outer Depth (01 May 98) 


106 


freq a -10 9.0 









Table G.17. LDV Data: Center Position, Inner Depth (01 May 98) 
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APPENDIX H. COMPUTED ROTOR EXIT PLANE DATA 


Pssst stance aera Te span [remaar | wever [apna Tp 1 psp J poor Ton Toe Tach coer | 


x 




























































































































































































































































































































































































































































































































































































































































































































































































































1 0.15536 -90 0.74998 | 0.7631 | 0.97177 | 0.9766 0.1576 | 0. 31483 
2 Rae ees 0.11679 0.75001 | 0.75742 | 0.96959 | 0.97232 | 0.11861 } 0.36273 
3 0.0002912 0.07242 0.01539 | 0.75006 | 0.75291 | 0.96852 | 0.96957 | 0.07359 | 0.39074 
4 0.0005413 0.05089 -0.03622 | 0.75015 | 0.75156 | 0.96604 | 0.96655 | 0.05177 | 0.42681 
5 0.0009042 0.08206 -0.13857 | 0.75025 — 0.75393 } 0.96229 | 0.96364 | 0.08365 | 0 46898 
6 0.0014302 0.12575 -0.30246 | 0.75033 | 0.75904 J 0.95789 | 0.96105 ] 0.12848 | 0.51423 
7 0.0021929 0.16482 -0.53341 | 0.75038 [ 0.76546 | 0.95345 | 0.95888 | 0.1688 | 0.55922 
8 0.0032975 0.197 -0.84602 | 0.75041 } 0.77211 7 0.94939 | 0.95715 | 0.20219 F 0.60153 
2 0.004897 0.22286 -1.25549 | 0.75043 | 0.77838 | 0.94589 | 0.95582 | 0.22915 | 0.63956 
0.0072085 0.24353 -1.76063 | 0.75047 | 0.78403 7 0.94299 } 0.95485 | 0.25079 | 0.67237 
0.0105442 0.26015 | 63.37898 | -2.29185 | 0.75055 | 0.78903 | 0.94064 | 0.95418 | 0.26823 | 0.70004 
0.0153442 0.27379 | 64.15511 | -2.68573 | 0.75069 | 0.7935 | 0.93873 | 0.95372 | 0.28258 | 0.72372 
0.0222257 0.28506 | 64.44872] -2.74014 | 0.75092 | 0.79749 | 0.93714 } 0.95339 | 0.29446 } 0.74446 
0.0320376 0.29419 { 64.18138 | -2.34168 | 0.75123 0.801 0.93586 | 0.95317 | 0.30411 | 0.76232 
0.0459186 0.30145 | 63.33211 ] -1.58692 | 0.75163 | 0.80403 | 0.93488 | 0.95305 | 0.31177 | 0.77703 
0.0653413 0.30734 | 61.9632 | -0.75162] 0.7521 [ 0.80669 | 0.93413 | 0.95302 | 0.31799 } 0.78923 
0.0921035 0.31257 | 60.2165 Ff -0.12544 | 0.75262 | 0.80922 | 0.93346 | 0.95301 | 0.32352 | 0.80075 
0.1282089 0.31765 | 58.27936 0.75322 | 0.81181 | 0.9328 | 0.95298 | 0.32889 | 0.81316 
0.1755558 0.32243 } 56.33136 0.75396 | 0.81448 | 0.93217 | 0.95296 | 0.33395 [| 0.82604 

0.2353854 0.32606 [| 54.50012 0.7549 § 0.81694 [| 0.9318 | 0.95306 | 0.33778 | 0.83614 

0.3075458 0.32716 | 52.84881 0.75608 | 0.81862 | 0.93202 | 0.95343 | 0.33888 | 0.83802 
0.3898437 0.32395 | 51.38651 0.75742 | 0.81875 | 0.93317 | 0.95416 ] 0.33535 | 0.82545 

0.4779231 0.31456 | 50.0649 0.75879 | 0.81647 } 0.93553 ] 0.95532 | 0.32522 | 0.79382 
0.5660036 0.29745 | 48.72149 0.75999 | 0.8113 0.9392 0.9569 | 0.30693 | 0.74314 
0.6483053 0.27243 | 46.93533 0.76097 | 0.80367 0.944 0.95884 | 0.28039 | 0.67989 

0.720471 : 0.24227 | 43.8221 0.76183 | 0.79531 | 0.94928 | 0.96102 § 0.24865 | 0.61535 

0.7803058 78.03058 0.21432 | 38.12332 0.76284 | 0.78886 | 0.9541 § 0.96329 | 0.21942 [| 0.56037 

0.8276587 82.76587 0.19864 | 29.54634] 4. 0.76409 | 0.78636 | 0.95766 | 0.96555 | 0.20298 | 0.51913 

0.8637684 86.37684 0.19822 | 20.30036 | 6.18799 | 0.76506 | 0.78721 } 0.95971 | 0.96757 | 0.20234 | 0.49087 
0.8905344 89.05344 | 85.68281 } 0.20316 | 12.46048] 6.78878 §| 0.7657 | 0.78897 | 0.96092 | 0.96917 | 0.20725 | 0.47079 

0.9099601 90.99601 | 87.49698 | 0.20452 | 5.74451 | 6.19076 | 0.76669 | 0.79028 —| 0.9624 | 0.97077 | 0.20848 | 0.44948 
0.9238429 92.38429 | 88.90458 | 0.20505 | -0.0977 | 4.92288 | 0.76759 | 0.79128 | 0.96389 | 0.9723 | 0.20885 | 0.42816 
0.9336562 4. 1 5661 | 93.36562 | 89.99778 | 0.20994 | 4.15661 | 3.7093 | 0.76783 | 0.79267 | 0.96469 | 0.97351 | 0.21375 | 0.41247 

0.9405392 -5.92609 | 94.05392 | 90.83712 | 0.21878 f -5.92609 | 2.9262 | 0.76765 | 0.79465 | 0.96468 | 0.97426 | 0.22274 | 0.40499 

0.94534 -6.02222 § 94.534 [| 91.46902] 0.22825 | -6.02222 | 2.53194 | 0.76737 § 0.79679 | 0.96408 | 0.9745 } 0.23246 | 0.40578 

0.9486767 -5.37274 1 94.86767 | 91.9342 | 0.23602 | -5.37274 | 2.37612 | 0.76711 | 0.79861 | 0.96326 | 0.9744 | 0.24048 | 0.4113 

0.950988 4.60695 | 95.0988 | 92.26965 | 0.24155 | 4.60695 | 2.3469 | 0.76691 | 0.79994 | 0.96257 | 0.97424 | 0.2462 | 0.41696 

0.9525878 -3.98802 | 95.25878 | 92.50778} 0.24531 | -3.98802 | 2.36877 | 0.76676 | 0.80086 | 0.96211 | 0.97414 | 0.25009 | 0.42058 
0.9536924 -3.56023 | 95.36924 | 92.67508] 0.24793 | -3.56023 | 2.39648 | 0.76667 | 0.80152 | 0.96183 | 0.97413 § 0.2528 | 0.42236 

0.9544551 -3.28473 | 95.44551 | 92.79187 | 0.24986 | -3.28473 | 2.41211 | 0.76661 | 0.80202 | 0.96165 | 0.97414 | 0.25479 | 0.4233 

0.9549811 -3.11026 | 95.49811 | 92.87312 } 0.2513 ]| -3.11026 ] 2.41676 | 0.76657 | 0.8024 § 0.96151 | 0.97414 | 0.25628 } 0.42415 
0.9553448 -2.99251 | 95.53448 | 92.92952 | 0.25235 § -2.99251 | 2.41856 | 0.76654 | 0.80268 | 0.96138 }] 0.97412 | 0.25737 | 0.42518 
0.9555948 -2.90457 | 95.55948 | 92.96933] 0.2531 | -2.90457 | 2.41968 | 0.76652 | 0.80288 | 0.96127 | 0.97408 | 0.25815 | 0.42626 
0.9557878 -2.84507 | 95.57878 | 92.99705 1 0.25359 | -2.84507 | 2.42149 | 0.76651 | 0.80301 | 0.9612 } 0.97406 | 0.25866 ]| 0.42699 
0.9558855 -2.79923 | 95.58855 | 93.01547 | 0.25389 | -2.79923 | 2.42504 | 0.7665 | 0.80309 | 0.96115 | 0.97404 | 0.25897 | 0.42746 
0.9559837 -2.74541 | 95.59837 | 93.04699 | 0.25421 | -2.74541 | 2.42885 | 0.76648 | 0.80317 | 0.96103 | 0.97395 | 0.25931 | 0.42909 
0.9565206 -2.56355 | 95.65206 § 93.13898] 0.25541 2.44154 | 0.76644 } 0.80349 | 0.96076 | 0.97381 | 0.26057 | 0.43216 
0.9577637 -2.03444 | 95.77637 } 93.37547] 0.25905 2.47181 | 0.76633 | 0.80449 | 0.96007 | 0.97349 § 0.26438 | 0.4399 
0.9608718 -0.71917 | 96.08718 | 93.93501 | 0.26733 2.61209 | 0.76611 | 0.80685 | 0.95843 | 0.97272 § 0.27307 | 0.45864 

0.9675 166 1.52727 | 96.75166 | 95.04134] 0.27939 2.77748 | 0.76564 } 0.81029 | 0.95622 | 0.97184 | 0.28571 | 0.48282 
0.977943 4.14735 | 97.7943 | 96.65221] 0.2877 2.71968 | 0.76495 § 0.81238 | 0.95474 | 0.9713 | 0.29444 | 0.49798 

0.9883689 6.6478 [| 98.83689 | 98.25707] 0.2814 2.43416 | 0.76429 | 0.80956 | 0.95547 | 0.97131 | 0.28789 § 0.48984 

0.9950138 8.71963 | 99.50138 | 99.32412 | 0.25256 2.07648 } 0.76376 | 0.7999 | 0.95932 | 0.97207 | 0.25786 } 0.45192 
0.9981218 9.99069 | 99.81218 | 99.80512 | 0.20055 1.81397 | 0.76315 | 0.78564 }| 0.96552 | 0.97357 | 0.2041 0.39682 
0.9993649 10.35494 | 99.93649 | 99.95979]} 0.1302 1.70004 } 0.76247 | 0.77183 | 0.97148 | 0.97487 | 0.13209 [| 0.35227 

0.9999018 9.8434 §99.99018 7 99.99555 }| 0.0561 1.68314 | 0.76201 | 0.76374 | 0.9741 | 097473 ] 0.05685 | 0.34099 








1 0 100 100 0 0 0.76135 { 0.76135 | 0.9791 0.9791 0 0.27883 








Table H.1. Rotor Exit Plan Data 
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APPENDIX I. SECONDARY PROGRAM 


The SECONDARY program takes the grid file and the solution file 
and creates a plane grid with a solution file. This is convenient for 


viewing a flat exit plan at a given location behind the leading edge of the 


blade. 


CE KEKE KAKKEREREKEREKEKEEKERKREKKEKEKEKEKEREKERERK KEK RKAEKKKKARRREKKRKKK KK KE 


Cc pxy.f reads rvc3d files & writes ascii files for plotxy 

e Mnit 1 = Anowe xyz fre 

C une os «= inputs q £ie 

c unit 6 = input error im-imax.or.jm-jmax.or.km-kmax.ne.0 

c unit 29 = output secondary.gq 

c unit 30 = output secondary.x 

CE RR KKK KKK KKK KEKE KEK KKK KEKE KK KEKKEEKKEKKEKKEKKEKEKKKKKEKEKKKKKKKKKK KK KKK KKK KKK KKK 


parameter (ni=350,nj=51,nk=72) 

rea lex<(ma, nj Mice y(n, Malo (nae, ra, ne) 

real qq({5,;nl,nj,nk);resda (500075) 

real u(ni,nj nk), v(m pe wie ak) 

real -xpp(2*nj7 nk) yop (2 hI mk) ;2Zpoi2 hank) 
real upp(2*nj7nk) 7 vent *ng4 me) wore (27h nk) 
real- xps (2 nig) yosts NI me) Zoster) 
real wps (2*ma7nk) vos (27° na nk) (wos (2 "ng, ak) 
real x2d(2*nj,nk) , y2d(2*nj,nk) ,z2d(2*nj,nk) ,q2d(5,2*nj,nk) 
Kealy xia 2 “1 Vy late clot mt Go oe 2s) 
Meal wlat2 ha 7 Vv lel2 ne ewe 2 ns) 

integer imd2,kmd2 


CE KKK KKK KEKE KKKEEKEKEKEKEEKEKEKEKKEKEKEKKEKEKKEEKRKEEKRKEKKEKKEKKEKKEKKKEKKEKKKEKKKKKKKKKKKKEKKK 


c read grid coordinates 
CER KKK KEK KEK KEKE KEKE KEEEKEKEKEKEKEKEREKEKKEKEKEKEKEKKEKKEKKKKKKKKKKKKEKKKKRKKKKKKKKKKEK 


read( L)am, jm, mM 
reads ( (GK) 2 = ame aii) a ke 1 em) 


a (iy Ais ga 40) a =i) Fs eT ee ee) 

2 (NCZ (aeeae, Ie) jaa en ee ea) eer 
CORE KKK KKK KEKE KEKE KEE KEKKEKEEKEKEK KKK KKKEKKEKKEHE KKK KKKKREKKKEKKEKRKKKEKEKKEKAEKKKEKKKKKEKSE 
Cc read restart file 


QE RR KKKKKEKEKKEKEKEKEEKEKKEKEKEEKEKEEKKEKEEKREKEKEKREKEKREKEKEEKEKEKERKEEKRKKKKKKKKKKKK K 


read(3) imax,  jJmax, kmax 
read(3)fsmach,alpha,re,time 


e 
icheck=iabs (im-imax) +iabs (jm-jmax) +iabs (km-kmax) 
1£(icheck.ne.0) then 
write(6)im, jm, km, imax, jmax, kmax 
stop 
endif 
c 
read (3),G0( (dat 13k) 4 2H mee om) ka oe oe 
c 
e additional residual data 
read(3)itl,iil1,phdeg,ga,om,nres,dum, dum, dum, dum 
read( 3s) weso Gr ye nmra=l1 nres) 71-1), 
e 
Drines *- Done reading Eort 3" 
e 
c non-dimensionalize the velocity field wrt Vref 


Til 


ImMaz=21mM7 2 

kmd2=km/2 
ul=qq(2,imd2, jm, kmd2)/qq(1,imd2, jm, kmd2) 
vl=qq(3,imd2, jm, kmd2)/qq(1,imd2,4jm,kmd2) 
wl=qq(4,imd2, jm, kmd2) /qq(1, imd2, jm, kmd2) 
vref=saqrt (ul**2+vl**2+wl**2Z) 


Dring = .! “Vrer = “,vrer 
CG 
Cc calculate the velocity field 
C 
GO71=2.- 1m 
go y= i am 
do k=1,km 


Wet 7. K=O a tek) do ae ik) 7 eee 
Via deka (35 ide <) ag (lea a yy vet 
Wl 35 e) OC (454.9), 6) 7 od (Ll, eee 
enddo 
enddo 
enddo 
print *,u (im, 7m, km) 7 v (im, jm, we a) 


CO I I IK KI II II KI III ITO III I OO Ok kk kk 
(ohhh hh hehe ee eee eee eee ee ee ee eee eS 2 ee 


C 
e 
Cc find the xmin and xmax grid points corresponding to the min and 
max 
C axial location of the blade 
e 
Sere iad alae 
xmax=0. 
kaa 
G 
PS tame=1t iL 
Rei t=im—1t 1 +1 
dois =I stant, ran & 
2£(x(i,1,k)..le.samin) then 
a 
Stix ( Uo ee) 
endif 
é 
1£(x(1,1,k).ge.xmax) then 
imax=1 
xmax=x(1,1,k) 
endif 
enddo 
Detie =) 'XMMN—=' sanin, sanax— cone 
C 
c find the grid no. at 130% chord 
c 
chord=abs (xmax-xmin) 
Pride) chord="chord 
Cilm=i-685147* echord+xmin 
Dette veh ys ern lie 
C 
ao. K=17500 
dena) 97m 


dowi= ll rmeaz— i 
ehillo—-x (1,3, Kk) 
1£(chlo.le.chlr) then 
C= sol Lak — x (1,50) 
Mppiag.. «) chi 


2 


On 


1000 


OO T 


Yop (a7 k)=y(2-1,9.h)> 
4 (oe 14 Kk) Sep ie yas a ie 1) 
Zpp (ay kK) =z (i=, 77k) 
Hix 114k) ene) (2 eZ aed 
Upp (47K) u(t a) 
(sc (i-dy 5 kK) Chl) 4 (ea ee a - 2, hood 
Vp (15) =v aa) 
FAM (=a kK) =Chil) = (tia kK) (2 ea Ke) 
wpp(j,k)=w(1i-1,3,k)+ 
4x ( ta 34) — Chl) Aw, hew a - bao k) ed 
Go =o 1000 
endif 
enddo 
do i=im,1imd2+1,-1 
chlo=x(1,j,k) 
1£ (Ehlo,. le emia) then 
Case (tl Soe le ala) 
SOS 1) Kk) emia 
yps(j,k)=y(i.j,k)+ 
#(ehir x(a eK) Ga ee ye a 
zps(j,k)=z(i,j,k)+ 
4{(cHilr=x (i759 ,K) ) 2 eee) 2 i) 
ups (j,k) =u(i,j,k)+ 
# (Chiyr—x(155,4)9* (U2 eel Ko 
VS (ay Kk) =v eg) > 
Hf (Ch 5 (a oa aa atl de ee ii) yea 
wps (j,k) =w(i,j,k)+ 
#{(Chir-=x(1, 35) (wii 9, kK) wilt, Kk) ya 
Go 20 1001 
endif 
enddo 
continue 
enddo 
enddo 


print *,"Finished interpolating" 


write plot3d 3d-file 


Weeder 2 Sad 2 Clea. Cl, k 
Write (30) 172a74 24,1 
Wwaedte( 3) 12d i L 

write (32)27Zd7e 

write (29) fsmach,alpha,re,time 
write(31)fsmach,alpha,re,time 


1—=0 
(deowj= 3m...) 
f= 1 
do k=1,km 
x2d(i,k)=xpp (j,k) 
y2d(i,k)=ypp(j,k) 
Zante k) =zZpp (a, Kk) 
GZ Gt ede) =e 
G26 274 2 e-Upp gi) 
q7at3 2 ki=vpp (i. 
q2a(4,1,k)=wpp (j,k) 
Gio in, ik) ke 
enddo 
enddo 


8. 


Go J=2; )m 
L=it+l 
ao k=i. km 
<2 Gi tek = xps (5 4) 
y2d(1i,k)=yps(j,k) 
ZI ada) =Zps ty.) 
Godt teak =i. 
Gz? 17k) =Upsita,) 
q2dt3, 1, )=vps (4.0) 
q2d(4,1,k)=wps(j,k) 
q2d(5,i,k)=1. 
enddo 
enddo 


measure=5.0 
ado 1=l,.42¢d 
do j=1, 32d 
E=sgretvyedCi iy ote ca ey a ae) 
if (measure.le.r)then 
d=y2d(i,j-1)-y2d(i,j) 
xd (a) =eh ia 
yld(i)=measure 
Zend) =2 20 eae 
# (y2d(1,j-1) -measure) *(z2d(1,3)-z2d(i,j-1))/d 
uld(i) =q2d(27175-1)- 
# (y2d(i,j3-1) -measure) * (q2d(2,1,j3)-q2d(2,i1,j-1))/d 
vild(1)=q2d(3,1-1,3)+ 
#(y2d(i,j-1) -measure) * (q2d(3,i1,j3)-q2d(3,i,j-1))/d 
wld(i)=q2d(4,1-1,j)+ 
#(y2d(i,j-1)-measure) * (q2d(4,1,j3)-q2d(4,i,j-1))/d 
qid(1,1)=1 
reply ea bllel (al) 
Calc ese ea) =v lech (.14) 
qld(4,i)=wl1d(i) 
emickGse.d ) =) 
Go to 20071 
endif 
enddo 
2004 continue 
enddo 
e 
Write (30) (x2d(1, 7) oad), j—1, 3:20) ,; 
# Gee cere), 1 2a a= 32a 
# ((z2a(i,j),i=1,i2a) ,j=1, 52d) 
Karta tee Oot (ep Cala etiee)) ye — lh lec) i ler Cleat lene 
write (32) ((xld(i),i=1,i2d),(yld(i),i=1,i2d), 
# Coc (a) = 1 ere) ) 
Wigdtee (Se chlch leila 1), 2c = lo) 
peint *,"FPinished” 


CREEK KKERKEKEKKKEREKEEKEKREEKEKEKEKKEKEKKEEEKEKEKKEEKKKKKEEKEKKKKKKEKKEKEKKKKKKEKEKEKE 


CEKKKKKEKEKEKKEKEEKKKEKEKEKEEKERKKEEKKKKEKKEEKKEKKEKEKEKEKKEKKEKEKEEKREREKRKEKKEKEKERKE 


CE KEK KK KEKE KK KEKE KEKE KEKE KKK KK KEKKKKEKEKKKEK KKK KKEKEKKAKKEKKKEKKKKKEKKKKKKKKKEKKKKKEKEK 


stop 
end 
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